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1.  INTRODUCTION 


Recent  developments  in  the  field  of  polymer  chemistry  are  being  used  increas¬ 
ingly  in  current  munition  design.  Increased  use  of  plastic  materials  is  due  to 
their  fabrication  flexibility,  high  strength-to-weignt  ratios,  and  low  cost. 
However,  in  munitions,  plastics  are  often  used  in  direct  contact  with  propel¬ 
lants  or  explosive  materials  To  ensure  cafely  in  handling  and  reliability 
after  .storage,  plastic  materials  must  be  compatible  with  propellant  mate¬ 
rials. 

Two  aspects  of  compatibility  must  be  considered.  The  plastic  must  not  react 
with  the  propellant,  creating  a  more  sensitive  or  degraded  product.  Similarly, 
the  propellant  must  not  attack  the  plastic,  resulting  in  degradation  of  the 
physical  properties  of  the  plastic.  Thermal  analysis  is  one  of  the  analytical 
test  techniques  useful  in  evaluating  the  stability  of  a  propellant  exposed  to  a 
plastic  material  during  the  operational  life  of  a  munition. 

The  test  program  described  in  this  report  was  performed  using  thermal 
analysis  techniques  to  evaluate  the  compatibility  of  propellants  IV1-26,  PYRO 
ano  NhCO  with  a  variety  ot  plastic  adhesives  and  loams,  me  data  generated 
give  a  statistical  measure  ol  propellant  stability  during  aging  in  contact  with 
the  plastics.  No  data  were  generated  regaining  the  effect  ol  the  propellant 
on  the  plastic  materials. 


II.  SUMMARY 


The  objective  of  this  program  was  to  evaluate  the  stability  of  several  pro¬ 
pellants  when  aged  in  contact  with  plastic  adhesives  and  foams.  Th  _?  data  are 
required  for  application  to  Navy  gun  propulsion  technology. 

There  are  several  tests  useful  in  determining  incompatibility  in  px  opellant- 
plastic  mixtures.  Rogers1  has  proposed  the  use  of  a  thermal  initiation  test  to 
determine  short-term  thermal  instability.  This  technique  is  a  modification  of 
the  older  Henkin  test2.  Long-term  aging  stability  has  been  predicted  prin¬ 
cipally  by  use  of  the  vacuum  stability  test3,  4  that  predicts  compatibility  by 
measuring  gas  evolution  of  the  propellant- plastic  system  when  aged  at  ele¬ 
vated  temperatures.  Variations  of  the  vacuum  stability  test  include  Taliani, 
modified  Taliani,  and  witness  tests.  These  variations  are  primarily  concerned 
with  differences  in  recording  and/or  analyzing  data. 

The  modified  Henkin  test  is  not  useful  for  predicting  long-term  aging  stability. 
The  vacuum  stability  test  will  indicate  short-term  aging  instability,  although 
the  results  must  be  interpreted  noting  that  typically  only  one  test  temperature 
is  used.  Several  limitations  should  be  noted  concerning  use  of  the  vacuum 
stability  test  to  determine  long-term  stability. 

First,  the  test  is  not  readily  adaptable  to  quantitative  analysis.  In  the  typical 
vacuum  stability  test,  as  performed  by  Picatinny  Arsenal  ,  a  net  increase  of 
five  milliliters  in  the  amount  of  gas  given  off  by  five  grams  of  propellant  (due 
to  exposure  to  the  plastic)  is  considered  evidence  of  excessive  reactivity.  The 
Naval  Ordnance  Laboratories  have  established  four  stability  classes  for  the 
same  test  (see  Table  I).  They  determined  each  class  by  the  volume  of  gas 
evolved  per  gram  propellant  after  48  hours  at  100°C  [(milliliters  of  gas 
measured  at  standard  temperature  and  pressure  (STF)]. 


Table  I.  Classification  of  Propellants  by  the  100°C  Vacuum  Stability  Te3t^ 


Class 

Milliliters  of  gas /gram/ 

48  hour  s/ ST P 

Quality 

I 

0-2 

Satisfactory 

Stability 

II 

2-6 

Borderline 

Stability 

III 

6-18 

Unstable 

IV 

Over  18 

Extremely  Unstable, 
Possibly  Hazardous 
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Evolution  of  additional  gas*  is  attributed  to  undesirable  interaction  between 
the  propellant  and  the  plastic,  or  undesirable  reactions  induced  by  the  poly- 
mer^.  The  gas  evolution  criterion,  to  be  quantitatively  useful,  must  be  cor¬ 
related  with  long-term  aging  tests  on  each  propellant-plastic  system.  This, 
of  course,  is  desirable  regardless  of  the  technique;  still,  it  is  not  completely 
practical  considering  the  large  number  of  propellants  and  plastics  available. 

To  be  of  practical  value,  a  test  technique  should  yield  data  that  will  allow  a 
kinetic  analysis  to  predict  aging  stability. 

The  second  limitation  of  the  vacuum  stability  test  is  that  the  polymers  are 
generally  tested  in  the  solid  phase  (thermoplastic  or  cured  thermoset).  In 
many  end  use  applications,  it  is  desirable  to  cure  the  plastic  in  contact  with 
the  propellant.  This  may  create  a  system  that  is  very  unstable  during  cure 
of  the  polymer,  but  is  compatible  after  the  polymer  i3  fully  cured.  This 
problem  can  be  overcome  by  testing  the  vacuum  stability  of  each  of  the  plas¬ 
tic  components  (resin,  curing  agent,  catalyst,  filler,  etc.  )  with  the  propel¬ 
lant  and  then  testing  the  resulting  cured  system.  Unfortunately,  in  some 
cases,  this  technique  does  not  duplicate  the  environment  experienced  by  the 
propellant,  and  it  also  markedly  increases  the  test  time. 

Thermal  analysis  test  techniques  such  as  Differential  Thermal  Analysis  (DTA), 
Differential  Scanning  Calorimetry  (DSC),  and  Effluent  Gas  Analysis  (EGA) 
can  be  used  to  determine  short-term  thermal  instability  and  to  predict  long¬ 
term  aging  stability.  Quantitative  predictions  of  long-term  aging  stability 
can  be  made  by  applying  one  of  several  kinetic  interpretations  to  the  data. 

Both  the  fully  cured  polymer  and  the  polymer  curing  against  the  propellant 
can  be  tested. 

Thermal  analysis  techniques  measure  enthalpic  effects  attendant  in  most 
chemical  or  physical  changes  that  occur  within  a  material  as  the  environmental 
temperature  is  increased  at  a  selected  rate.  During  a  thermal  analysis 
test,  the  temperature  of  the  sample  is  compared  with  the  temperature  of  a 
reference  material  while  both 'are  heated  at  a  controlled  rate.  The  results 
are  normally  plotted  as  temperature  difference  versus  reference  termperature. 
DSC  combines  the  DTA  method  of  temperature  comparison  with  closed-loop 
heat  compensation.  The  measurement  plotted  is  the  electrical  energy  in 
calories  per  second  required  to  maintain  the  sample  and  reference  at  the 
same  temperature.  The  results  plotted  show  endothermic  and  exothermic 
peaks  whose  areas  are  directly  proportional  to  the  energy  absorbed  or 
liberated  by  the  sample.  EGA  measures  the  relative  quantity  of  gas  evolved 
by  the  sample  as  the  environmental  temperature  is  increased  at  a  selected 
rate. 


’^Additional  gas  refers  to  the  quantity  of  gas  evolved  in  excess  of  that  evolved 
by  propellant  and  plastic  when  separately  exposed  to  48  hours  at  100°C. 
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In  using  thermal  analysis  to  evaluate  the  compatibility  of  propellant- plastic 
combinations,  the  assumption  is  made  that  either  sample  temperature  or  gas 
evolution  is  maximum  at  the  peak  reaction  rate.  This  assumption  is  well 
founded  for  propellants  and  explosives  as  indicated  by  Cook^. 
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III.  DISCUSSION 


Compatibility  of  the  test  plastics  with  various  gun  propellants  was  examined 
by  thermal  analysis  techniques,  primarily  Differential  Scanning  Calorimetry. 

Compatibility  testing  by  thermal  analysis  examines  the  kinetics  of  the  thermal 
decomposition  reaction.  The  compatibility  of  a  material  is  determined  by  the 
effect  exposure  to  the  material  has  on  the  decomposition  reaction  kinetics  of 
the  explosive  or  propellant. 

It  should  be  noted  that  these  tests  make  no  inferences  as  to  the  chemistry  of 
the  decomposition  reaction  or  to  the  effect  of  the  propellants  on  the  test  plastic 

Activation  energy  and  frequency  factor  of  the  decomposition  reaction  were  de¬ 
termined  by  the  variable  heating  rate  method  as  developed  by  Kissingerb  and 
modified  by  Pakulak7. 

The  reaction  kinetics  as  described  by  Kissinger  are  based  upon  first  order 
reaction  rates.  Assuming  a  first  order  reaction  rate,  an  explosive  will  de¬ 
grade  as  follows: 

|if)T=  *T  (1-X)  (1) 

where  x  =  the  fraction  of  the  explosive  degraded 
t  =  time 

k  =  the  reaction  rate  constant  at  temperature  (T). 

The  reaction  rate  constant  is  determined  from  the  Arrhenius  Rate  Equation 

k(T)  =  A  exp  (-E/RT)  (2) 

where  A  =  frequency  factor 
E  =  activation  energy 

R  =  gas  constant  (1.  987  calories /mole°K) 

T  =  temperature  in  degrees  Kelvin. 

When  the  temperature  is  changing,  as  in  a  Differential  Scanning  Calorimeter 
run,  the  reaction  rate  becomes 
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dx  _ 
dt 


+ 


dx  _dT 
dT  t  dt 


(3) 


By  determining  X  and  T  at  the  same  instant,  the  time  coordinate 
becomes  zero.  This  reduces  the  rate  equation  to 


t 


~  =A(l-x)exp  (-E/RT) 


(4) 


When  the  reaction  is  occurring  at  its  maximum 
spect  to  time  becomes  zero.  Solving  for  I  ^  l 


rate, 

dx 
dt 


its  derivative  with  re- 
the  equation  becomes: 


_d_ 

dt 


dx 

dx 

(  E 

dt 

dt 

\  RT2 

=  A  exp  (-E/RT) 


(5) 


The  temperature  at  which  a  reaction  occurs  at  its  maximum  rate,  dx/dt,  can 
be  measured  as  the  peak  exotherm  temperature  of  a  DSC  thermogram.  There¬ 
fore,  the  maximum  reaction  rate  is  defined  as 


A  exp  (-E/RTm) 


dT 

dt 


16) 


This  equation  reduces  to 


(7) 


where  cp  is  the  change  in  temperature  with  respect  to  time  dT/dt  or  the  sample 
heating  rate. 

The  peak  exotherms  and  heating  rates  from  a  series  of  DSC  tests  can  be  used 
to  solve  for  the  activation  energy  by  Equation  (7).  Using  this  value  for  E, 
the  frequency  fact.  A,  can  be  determined  by  Equation  (6).  With  the  values 
of  E  and  A  determined,  the  rate  constant,  K,  can  be  determined  by  using  the 
Arrhenius  Rate  Equation,  (2). 


where  Tm  is  the  peak  temperature  of  a  DSC  trace, 
the  relationship 
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There  are  a  number  of  assumptions  necessary  in  making  a  prediction  of 
kinetic  data  from  thermal  analysis  data. 

1.  It  is  assumed  that  the  decomposition  by-products  are  a  combination 
of  gaseous  and  solid  nonexplosive  materials,  free  to  escape  without 
pressure  buildup.  An  isothermal  increase  in  pressure  due  to 
trapped  gases  would  affect  the  decomposition  reaction.  It  should 
be  noted  that  the  vacuum  stability  test  is  dependent  on  the  formation 
of  gaseous  by-products  for  a  prediction  of '  stability.  Thermal  de¬ 
composition  of  explosives  has  been  discussed  fry  several  authors, 
including  Gray  and  Bonomo8  and  Dacons,  et  alf. 

It  is  assumed  that  thermal  decomposition  values  are  comparable 
for  the  explosive  in  the  melt,  solid  and  vapor  phases.  This  as¬ 
sumption  is  incorrect,  but  the  difference  in  the  decomposition  rate 
for  the  two  phases  is  not  generally  known.  Various  authors 
claim  that  the  decomposition  rates  may  be  as  much  as  50  to  100 
times  faster  in  the  melt  phases  than  inthe  solid  phaoes.  For 
example.  Gray  and  Bonomo^  have  stated  that  production  of  TNT 
decomposition  by-products  is  dependent  on  temperature,  the 
greatest  decomposition  occurring  between  150°  and  160°F  (the 
highest  temperature  studied  in  their  work).  The  difference  is, 
of  course,  due  to  the  greater  mobility  of  the  atoms  involved  in 
the  decomposition  when  the  explosive  is  in  the  liquid  phase  (aboye 
its  melting  temperature).  Several  authors,  including  Farmer10 
and  Rosen1 1,  have  discussed  this  area. 

The  techniques  used  for  deriving  the  reaction  kinetics  from  thermal  analysis 
data  are  subject  to  much  discussion.  Reed,  et  al,  have  argued  that  Kiss¬ 
inger's  attempt  to  prove  that  the  maximum  reaction  rate  occurs  at  the  peak 
of  the  DTA  curve  is  incorrect.  They  note,  additionally,  that  regardless  of 
the  proof,  large  errors  in  the  calculated  value  of  E  can  occur  in  systems  when 
the  maximum  reaction  rates  are  within  1°C  (experimental  error)  of  the  peak 
temperature.  Masting  in  a  theoretical  study,  has  compared  Kissinger' s 
method  with  the  modification  used  by  Pakulak.  In  Kissinger's  original  method, 
the  actual  sample  temperature  at  the  peak  differential  temperature  is  used  in 
Equation  (7).  The  modified  formula  involves  the  use  of  the  reference  temper¬ 
ature  at  the  maximum  reaction  rate.  Mastin1^  concludes  that  the  Kissinger 
method  will  give  reasonably  good  kinetic  rate  parameters  in  a  simple  decom¬ 
position  reaction.  He  notes  that  the  greatest  error  occurs  at  fast  heating 
rates  when  it  is  assumed  that  the  rate  of  change  of  the  sample  temperature  at 
the  maximum  reaction  rate  is  equal  to  the  apparatus  heating  rate. 

With  a  temperature- compensating  design,  as  in  the  Perkin- Elmer  DSC  ap¬ 
paratus,  this  problem  is  diminished;  however,  thermal  gradients  within  the 
sample  are  still  a  source  of  possible  error.  For  this  reason,  sample  sizes 
are  kept  to  a  minimum. 
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Given  theprocessing  conditions  (the  time  at  which  the  system  will  be  held  at 
the  processing  temperature),  an  arbitrary  criterion  can  be  established  for 
short-term  thermal  stability.  For  example,  such  a  criterion  might  require 
that  the  presence  of  test  material  with  the  explosive  increase  the  amount  of 
explosive  reacted  by  less  than  one-tenth  of  one  percent. 

It  is  somewhat  more  difficult  to  establish  a  criterion  for  long-term  aging 
stability  because  of  the  unknowns  in  storage  time  and  temperature.  In  this 
case,  a  safe  criterion  would  require  that  there  be  no  significant  increase  in 
the  decomposition  rate  due  to  *'  presence  of  the  test  material.  This  require¬ 
ment  dictates  that  the  activati .  <  ••no  jy  and  frequency  factor  for  the  propellant 
only  be  compared  with  the  same  tie  parameters  for  the  propellant  exposed 
to  the  plastic.  Using  the  Arrhemu  Rate  Equation,  predictions  may  be  made 
as  to  the  life  of  test  and  control  samples.  If  these  predictions  are  statistically 
similar,  the  system  is  compatible. 
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IV. 


TECHNICAL  APPROACH 


All 'testing  was  performed  with  a  Perkin-Elmer  model  DSC-1B  Differential 
Scan*ing  Calorimeter.  Sample  preparation  was  as  follows: 

Samples  weighing  approximately  5  milligrams  were  cut  from  each  of 
the  propellants  and  run  as  controls.  Then  similar  samples  were  run 
in  contact  with  each  of  the  various  test  plastics.  With  the  liquid  plas¬ 
tic  samples,  the  propellant  sample  was  coated  and  allowed  to  dry. 

With  the  solid  and  foam  samples,  the  plastic  was  either  cut  into  thin 
Aices  or  ground  to  a  powder  and  mixed  with  the  propellant  sample. 

The  samples  were  placed  into  aluminum  sample  pans  with  an  aluminum 
c«ver  crimped  over  the  top  of  the  sample.  The  sample  pans  were  placed 
in  the  DSC  sample  cells  and  covered  with  an  aluminum  dome.  The  temp¬ 
erature  at  which  the  decomposition  rate  reached  its  maximum,  the  peak 
eocotherm  on  the  DSC  thermogram  (Figure  1),  was  then  recorded. 

Due  to  niomlinearities  in  thermocouple  output,  the  indicated  temperature  and 
the  actual \temperature  of  a  DSC  cell  may  not  be  the  same  across  the  entire 
temperature  range  of  the  instrument.  Because  of  the  need  for  precise  temper¬ 
ature  accuracy,  it  was  necessary  to  calibrate  the  DSC  apparatus.  This  was 
done  by  determining  the  indicated  melt  points  of  known  standards  across  the 
entire  temperature  range  and  at  each  of  the  test  heating  rates.  This  operation 
was  performed  twice,  first  prior  to  the  initial  series  of  testing  and  again 
prior  to  a  series  of  requested  i  v  tests  on  some  of  the  data.  The  calibration 
data  is  listed  iff  Table  II.  The  temperature  correction  data  at  each  heating 
rate  was  found  to  be  a  second  order  curve.  The  coordinates  of  these  curves 
were  calculated  to  a  least  squares  fit  by  computer  and  placed  in  a  FORTRAN 
program  that  corrected  the  temperature  error  of  the  DSC  runs  and  printed 
this  data  into  data  files.  This  program  is  listed  in  Appendix  A  as  Corect— 
Temperature  Correction  Program. 

The  corrected  data, from  each  propellant  and  propellant-plastic  combination 
was  then  run  through  a  second  FORTRAN  program  to  calculate  the  activation 
energy,  frequency  factor,  and  reaction  rate  constant  at  any  desired  temper¬ 
ature,  This  program  will  also  predict  the  extent  of  the  decomposition  reac¬ 
tion  over  any  desired  period  of  time  at  any  desired  temperature.  This  is  an 
extrapolation  of  data  over  several  magnitudes  of  time  via  an  Arrhenius  Rate 
Equation,  which  inust  be  remembered  when  the  predicted  service  lifes  are 
examined.  This  program  is  listed  in  Appendix  A  as  Energy — Activation 
Energy  Program, 

One  criterion  for  judging  a  plastic  to  be  compatible  with  a  propellant  was  that 
there  be  no  statistically  significant  difference  between  the  control  and  test 
sets  of  data.  Ail  statistical  analysis  was  performed  at  95  and  99  percent 
confidence  levels. 
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MATERIAL 


TRUE 

MELT  POINT 
CK) 


INDICATED  MELT  POINT  CK)  AT 
HEATING  RATE  CC/MINUTE)  OF 


INITIAL  TESTING  -  JULY  1972 


IODINE  METAL 
INDIUM  METAL 
AMMONIUM  NITRATE 
EUTECTIC  SOLDER 
TIN  METAL 
LEAD 


394.2 

388.8 

386.0 

Em 

383.3 

383.3 

438.0 

*33.0 

430.8 

428.9 

428.5 

454.2 

448.0 

444.6 

442.6 

442.1 

463.0 

457.8 

455.2 

IfiSIS 

453.0 

452.6 

517.9 

511.9 

509.2 

508.0 

507.0 

506.7 

609.3 

604.4 

601.8 

600.2 

599.7 

599.5 

RETESTING  -  FEBRUARY  1973 


NAPTHALENE 

353.4 

354.2 

BENZOIC  ACID 

395.6 

399.0 

INDIUM  METAL 

429.8 

436.0 

2-CHL0R0ANTHRAQUIN0NE 

482.1 

489.1 

TIN  METAL 

505.1 

512.8 

ANTHRAQUINONE 

557.4 

565.7 

LEAD 

600.7 

605.4 

343.4 

343.5 

389.5 

389.4 

425.2 

425.5 

479.4 

479.8 

The  £irst  test  was  an  analysis  to  determine  if  the  plot  of  1/T  versus  log 
/ T“,  from  which  the  activation  energy  was  calculated,  was  a  fit  to  a 
straight  line.  This  was  done  by  the  FORTRAN  program  listed  in  Appendix  A 
as  PROG1 —Statistical  Analysis  -  Program  1. 

The  second  program,  listed  as  PROG2— Statistical  Analysis  -  Program  2, 
was  used  to  compare  the  control  and  test  data  for  each  propellant-plastic 
combination.  This  program  performed  three  tests.  First,  it  determined 
if  the  variances  in  each  data  set  v  ere  comparable.  Second,  it  determined 
if  the  plots  of  1/T  versus  log  cp / rf2  for  both  the  control  and  test  sets  of  data 
were  parallel.  Third,  it  determined  if  the  two  plots  were  coincidental. 

Each  of  these  statistical  programs  must  be  satisfied  in  order  that  meaning¬ 
ful  statements  may  be  made  regarding  the  compatibility  of  a  given  propellant- 
plastic  combination. 

The  plot  of  1/T  versus  log  cp  /T2  measures  the  order  of  the  reaction.  If  this 
plot  fails  to  fit  a  straight  line,  the  assumption  that  the  reaction  kinetics  are 
first  order  is  incorrect.  None  of  the  propellants  evaluated  failed  this  test. 

Homogeneity  of  the  sample  variances  is  determined  to  test  the  assumption 
that  the  sample  data  comes  from  the  same  population.  The  statistical  tech¬ 
nique  used  to  compare  the  parallelism  and  contingency  of  two  sets  of  data 
(propellant  only  versus  propellant-plastic  combination)  is  dependent  upon 
whether  or  not  the  variances  are  equal. 

The  condition  of  two  sets  of  data  being  parallel  determines  whether  or  not 
exposure  of  the  propellant  to  the  plastic  has  significantly  changed  the  acti¬ 
vation  energy  of  the  propellant  decomposition  reaction..  Under  normal 
circumstances,  the  plastic  will  not  change  the  activation  energy  of  the  pro¬ 
pellant  regardless  of  compatibility.  When  it  occurs,  a  change  in  activation 
energy  may  indicate  an  incompatible  propellant -plastic  combination. 

The  condition  of  two  sets  of  data  being  contingent  determines  whether  expo¬ 
sure  of  the  propellant  to  the  plastic  has  significantly  changed  the  frequency 
factor  of  the  propellant  decomposition  reaction.  Under  normal  circumstances, 
for  a  propellant -pi  a  Stic  combination  that  is  not  compatible  the  activation 
energy  will  not  be  changed  (the  lines  will  be  parallel),  but  the  frequency  factor 
will  be  significantly  increased  (the  lines  will  not  be  contingent).  In  the  event 
of  the  frequency  factor  decreasing,  it  may  be  said  that  there  is  interaction 
between  the  plastic  and  propellant.,  but  the  result  is  not  such  that  a  degraded 
or  hazardous  propellant  will  result. 

The  results  of  the  testing  list  any  plastic-propellant  combination  that  passed 
all  of  the  tests  as  compatible.  If  a  plastic -propellant  combination  failed  at 
one  confidence  level  and  passed  at  the  other,  or  failed  any  of  the  statistical 
tests  at  both  confidence  levels  it  was  analyzed  further  by  reviewing  the  shelf 
rediction  obtained  from  the  Arrhenius  Rate  Equation.  Only  in  the  case 
wn  a  plastic -propellant  combination  failed  the  statistical  tests  and  indicated 

a  lat,  change  in  shelf  life  was  the  system  considered  incompatible. 
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V.  RESULTS  AND  CONCLUSIONS 


The  results  of  the  DSC  tests  on  each  propellant  and  propellant -plastic  combi¬ 
nation  are  included  in  Appendix  B.  The  temperatures  listed  are  the  corrected 
values.  Included  are  the  predictions  for  the  percent  of  propellant  unreacted 
when  stored  for  1,  2,  5,  10,  and  20  years  at  25°,  50°,  75",  and  100°C  (298°, 
323°,  348°,  and  373°K).  Also  attached  are  the  plots  of  1/T  versus  cp  /t^  for 
control  and  propellant-plastic  combination. 

Table  III  lists  the  results  of  the  statistical  analysis  of  the  test  data. 

The  data  fits  on  NACO  and  PYRO  control  and  plastic -propellant  combinations 
were  very  good.  The  data  for  M-26  was  much  more  scattered.  For  this' 
reason,  the  limits  for  statistical  noninteraction  for  M-26  tended  to  be  quite 
broad,  whereas  the  limits  for  the  other  propellants  were  more  exacting. 

This  accounts  for  the  fact  that  some  of  the  M-26  plastic  combinations  did  not 
statiscally  interact,  while  their  plots  seemed  to  be  more  divergent  than  some 
of  the  other  propellant-plastic  combinations  where  statistical  interaction 
was  noted. 

Several  combinations  are  noted  where  statistical  propellant -plastic  inter¬ 
actions  are  noted.  Reviewing  the  shelf-life  predictions  for  each  of  these 
combination  yields  two  examples  of  the  plastic  having  a  significant  effect 
on  the  propellant.  In  each  instance,  the  data  would  appear  to  suggest  greater 
stability  for  the  propellant  due  to  exposure  to  the  plastic.  It  is  concluded 
that  regardless  of  the  effect  predicted  by  extrapolation  with  the  Arrenius 
Rate  Equation  the  systems  have  demonstrated  a  large  chemical  interaction 
that  could  be  deleterious  to  device  performance.  For  this  reason  M-26  with 
EC  1099  and  PYRO  with  Ultrabond  76-125  are  considered  incompatible.  In 
all  other  cases,  the  combinations  are  judged  compatible. 

Compatibility  data  was  requested  on  Celcon  1011  with  M-26  propellant.  Con¬ 
siderable  difficulty  was  experienced  in  attempting  to  achieve  good  DSC  data 
due  to  melting  of  the  Celcon  prior  to  the  decomposition  of  M-26.  As  an  alter¬ 
native  approach,  the  Effluent  Gas  Analysis  (EGA)  accessory  to  the  DSC  was 
used  in  an  attempt  to  determine  the  M-26  decomposition  temperatures.  EGA 
measures  the  relative  quantities  of  gas  output  of  a  sample  during  a  DSC  run. 
Assuming  that  the  peak  decomposition  temperature  and  the  temperature  at 
which  peak  gas  evolution  occurs  are  the  same,  EGA  data  should  coincide  with 
DSC  data. 
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During  the  initial  EGA  testing  with  M-26  and  Celcon,  two  problems  were  en¬ 
countered.  First,  a  short  time  lag  exists  between  the  evolution  of  the  gas 
output  of  the  sample  and  its  detection  by  the  EGA  detector.  This  is  due  to 
time  required  for  the  gas  to  pass  from  the  sample  cell  through  a  line  to  the 
EGA  detector.  Secondly,  the  EGA  is  presently  set  up  to  operate  with  nitrogen 
as  a  carrier  gas.  The  use  of  a  nitrogen  atmosphere  rather  than  air  causes 
some  changes  in  the  decomposition  characteristics  of  the  M-26  propellant. 
These  problems  could  be  eliminated  for  future  testing  by  careful  calibration 
of  the  EGA  time  lag  and  conversion  to  air  as  a  carrier  gas. 


Table  III.  Statistical  Analysis  Results 


TEST  SAMPLE 

ACTIVATION 

ENERGY 

(KCAL/MOLEI 

— 

FREQUENCY 

FACTOR 

PROG  ONE 

VARIANCES 

PARALLEL 

COINCIDENT 

RATING 

M-26  CONTROL 

49.3 

0.75  X  1021 

P 

... 

... 

— 

9  <2 

p 

P 

P 

P 

COMPATIBLE 

M-26/GUN  PLUG  FOAM 

53.9 

o.n  x  io^ 

M -2  6/GUN  PLUG  FOAM  (HUMIDITY) 

54.6 

0.27  X  1024 

P 

P 

P 

P 

COMPATIBLE 

M-26/L0CTITE  A 

49.8 

0.15  X  1622 

P 

P 

P 

P 

COMPATIBLE 

M-26/L0CTITE  N 

47.4 

0.10  X  1021 

P 

P 

P 

p 

COMPATIBLE 

M-26/EC  1099 

59.6 

0.81  X  102h 

P 

P 

P 

F 

INCOMPATIBl  F 

M-26/UNFILLE0  CARTRIDGE  CASE 

49.7 

o.io  x  io22 

P 

P 

P 

P 

COMPATIBLE 

M-26/FILLED  CARTRIDGE  CASE 

47.7 

0.12  X  1021 

P 

P 

P 

P 

COMPATIBLE 

PYRO  CONTROL 

44.4 

0.19  X  1019 

P 

... 

... 

... 

pyro/unfilleo  cartridge  case 

45.8 

0.93  X  1019 

P 

P 

M 

M 

COMPATIBLE 

pyro/filleo  cartridge  case 

44.6 

0.26  X  1019 

P 

P 

P 

P 

COMPATIBLE 

PYRO/ULTRABONO  76-125 

52.6 

0.20  X  1023 

P 

P 

M 

M 

INCOMPATIBLE 

P 

P 

F 

COMPATIBLE 

PYRO/EC  1099 

46.6 

0.2j  X  IO*” 

NACO  CONTROL 

47.4 

0,64  X  1020 

P 

... 

— 

... 

NACO/EC  1099 

43.7 

0.13  X  10l 3 

P 

P 

M 

M 

COMPATIBLE 

NACO/ULTRABONO  76-125 

44.6 

0.36  x  1019 

P 

P 

P 

P 

COMPATIBLE 
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-16- 


CORECT 

FORTRAN  -DSC  TEMPERATURE  CORRECTION 


210  READ<1,50)LN0,P<N),TEMP 
220  50  FORMA T( 15,8010.0) 

230  WRITEC9,60)P(N) 

240  60  FORMAT(//’PHI:* ,F7 .3/ ‘MEASURED  CALCULATED*) 

250  DO  70  1=1,8 

260  IF(PARAM(I)-P(N))70,80,70 

270  70  CONTINUE 

280  STOP  PHI  NOT  IN  LIST 

290  80  DO  90  J=l,7 

300  IF(TEMP(J))100,100,90 

310  90  CONTINUE 

320  J=8 

325  100  J=J-1 

330  I NDEX( N) = J 

340  DO  120  K= 1 , J 

350  T( N,K) =A( 1,1) 

360  DO  110  L= 1 ,2 

370  110  T(N,K)=T(N,K)*TEMP<K)+A(L+1 ,1) 

380  120  WRITE(9, l 30) TEMP <K) ,T(N,K) 

390  130  F0RMAT<Fs.2,F12.2) 

400  G0T040 
410  140  LNC= 1 
420  WRITE<2, 1 50)LN0, N 
430  150  FORMA T( I4,1H, ,14) 

440  DO  160  1=1, N 
450  LN0=LN0+1 
460  J= I NDEXC I ) 

470  IF0R(8) =MFOR< J) 

480  160  WRITE (2, IFOR )LN0,P(I),CT(I,K),K=1,J) 

490  STOP  SAVE  YOUR  FILE 
500  END 
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ENERGY 

SuHTHAN  -  ACTIVATION  ENERGY  PROGRAM 


10  COMMON  X(50,2),;';.5Q)  ,A(2,?)  ,IJ(7)  ,CC50)  ,T(9)  ,4K(9)  ,Z(y  ) 

20  +,  0(9)  ,B(2>  ,D£G(9)  .NAPHG)  ,AK  1  (S)  ,EX(  50  ,WYE(50> 

30  DOUBLE  PRECISION  X , Y , A ,C , AK , 0 , H , YOU , DVEE , EO ,E , R , D , Z EE 
40  +,TEE,EFF,F.Sl ,S2,S3,EL,DF,P 
50  COMPLEX  RATIO 
60  RATIO- * 1NFI NIIE' 

70  WR I  T E  (  9 , 1 0 ) 

80  10  FORMA T( /* WANT  DATA  KILE  DESCRIPTION?  (YES  OR  NO)’) 

90  READ (9, 20)1  WANT 

100  20  FORMA  T ( SAC ) 

tlC  If ( I  WANT . EQ ,2HN0) OOTOAU 

120  WRITE 0,30) 

130  30  FORMAT!/ ’ThIS  PROGRAM  KITS  A  1-ST  DEGREE  POLYNOMIAL  TO  A  BIVAril’ 
i40  +  ’ATE’/’SET  OF  DATA,  AND  PRINTS  A  SUMMARY  OF  THE  FIT. ’/’THE’ 

150  +  •  DATA  IS  READ  FROM  A  FILE  WHICH  MUST  HE  PREPARED  PRIOR  ’ 

ISC  +  'TO  EXECUTION. V  THE  DATA  FILE  SHOULD  BE  SET  UP  AS  FOLLOWS:’// 

170  +  ’LINK  I:  LINE  N0,M*/’L1NE  2:  LINE  NO ,  V 1  ,  U  1  ,  U2  ,  U  3  ,  U  4  ,  U5 , 116  ,  U  7  *  / 

130  «  'LINE  3:  LINE  NO  ,  V2  ,  U  I  ,  U2  ,  IJ3  ,  U  4  ,  US  ,  US  ,  U7  ’  /  ’ . ETC’ 

190  +  ’  UNTIL  ALL  V-U  DATA  SETS  HAVE  BEEN  ENTERED V WHEHt  M  IS  ' 

200  +  ’THE  NUMBER  OK  V, U  DATA  SETS.’/’VN  IS  A  HEATING  RATE  IN  ’ 

210  +  'DEGREES  C.  PER  MIN.  ’  /  ’  <•  N  0  EACH  UN  IS  A  TEMPERATURE  IN  JturtfctJ  ’ 

220  +  ’KELVIN’/’  OBSERVED  AT  THE  HEATING  RATE  VN.V 

23 C  -i  ’UP  TO  7  TEMPERATURES  MAY  BE  USED  FOR  1  HEATING  RATE.’/ 

24U  +  ‘EACH  LINt  OF  THE  DATA  FILE  MUST  BEGIN  ’ 

250  +  ’WITH  A  LINE  NUMBER ’ / ’ FOLLOWED  BY  a  COMMA. ’//’ EACH  ENTRY  ’ 

260  -r  ’EXCEPT  THE  LAST  0»  ANY  LINE  SHOULD  BE  FOLLOWED  BY  A  COMMA.’// 

270  +  ’FOLLOW  ThE  ABOVE  P  AT  A  WITH  LINES  CONl A 1 N I N 0 ’ / ’ ( 1 N  ADDITION’ 

230  +  ’  TO  THE  LINE  NUMBER 

290  +  ’  A)  THE  NUMBER  OF  STORAGE  TEMPERATURES.  (ZERO  IF  NONE.)’/ 

300  +  ’  *  IGNORE  STtPS  B.O.D.AND  E  IF  THE  NUMBER  ABOVE  IS  C.»*’/ 

310  +  '  B)  THE  PARTICULAR  STORAGE  TEMPERATURES  ( DEGRFES  KELVIN)’/ 

320  +  ’  C)  THE  NUMBER  OF  VALUES  FOLLOWING  (ZERO  IF  NONE.)'/ 

330  +  '  D)  EITHER  Cl/CO  OH  TIME,  WHICHEVER  IS  TO  BE  SOLVED  FOR.’/ 

340  +  ’  E)  PARAMETER  VALUES  FOR  TIME  (YEARS),  OR  Cl/CO  (PERCENT)’ 

350  +  /’  -DEPENDING  UPON  WHAT  WORD  WAS  TYPED  FOR  PARI  D,*//) 

360  STOP 

370  40  WR I TE( 9 , 50  ) 

3gC  50  FORMAT! /’NAME  OF  INPUT  DATA  FILE’) 

390  REA D ( 9 , 20 ) NAM 
400  CALL  DEFI NEC  t , NAM ) 

410  WRITEC9.70) 

420  70  FORMAT!/ ’NAME  OF  OUTPUT  DA T«  FILE’) 

430  READ! 9 ,20) NAM 
440  CALL  DEFI NE (2 , NAM ) 

450  LNO-0 

460  READ!  I , 80 ) LNS , M 
470  80  FORMAT! 2 16) 

480  S3  :0  .  DO 
490  WRITE (9, 90) 

500  90  FORMAT! //’INPUT  DATA’/’  PHI  TEMP ’ 9 X ’ I / T ’ I  1 X ’ PH  I / 1 SOU ARE  *  ) 
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ENERGY 

FORTRAN  -  ACTIVATION  ENERGY  PROGRAM 


10  COMMON  X(50,2),Y(50),A(2,2),U(7),C(50),T(9),AK(9),Z(9) 
20  4-,  Q<9) ,B(2) ,DEG(9) , NAM(6) ,AKI (9) ,EX(50> ,WYE(50) 

30  DOUBLE  PRECISION  X , Y  ,A,C ,AK ,Q ,B, YOU.DVEE, EO ,E,R, D ,ZEE 
40  +,TEE,EFF,F,S1,S2,S3,EL,DF,P 
50  COMPLEX  RATIO 
60  RATIO: 'INFINITE' 

70  WRITE(9, 10) 

80  10  FORMAT(/'WANT 
90  READ(9,20)IWANT 
100  20  F0RMAT(6A2) 


DATA  FILE  DESCRIPTION?  (YES  OR  NO)’) 


110  IF( I WANT.EQ.2HN0) G0T040 
120  WRITE(9,30) 

130  30  FORMAT(/*THIS  PROGRAM  FITS  A  1  -ST  DEGREE  POLYNOMIAL  TO  A  BIVARI' 
140  +  *ATE'/*SET  OF  DATA,  AND  PRINTS  A  SUMMARY  OF  THE  FIT. '/’THE' 

150  +  *  DATA  IS  READ  FROM  A  FILE  WHICH  MUST  BE  PREPARED  PRIOR  * 

ISO  +  'TO  EXECUTION. '/'THE  DATA  FILE  SHOULD  BE  SET  UP  AS  FOLLOWS:'// 

170  +  'LINE  1:  LINE  NO,M*/'LINE  2:  LINE  NO, VI ,U! ,U2,U3,U4,U5,U6,U7*/ 

180  +  'LINE  3:  LINE  NO, V2,U 1 ,U2,U3 ,U4,U5 ,U6 ,U7 '/’ . ETC' 

190  +  '  UNTIL  ALL  V-U  DATA  SETS  HAVE  BEEN  ENTERED V ' WHERE  M  IS  ' 

200  +  'THE  NUMBER  OF  V,U  DATA  SETS,'/*VN  IS  A  HEATING  RATE  IN  ' 

210  +'DEGREES  C.  PER  MIN., '/'AND  EACH  UN  IS  A  TEMPERATURE  IN  DEGREES  ' 

220  +  'KELVIN'/'  OBSERVED  AT  THE  HEATING  RATE  VN. V 

230  4-’ UP  TO  7  TEMPERATURES  MAY  BE  USED  FOR  1  HEATING  RATE.’/ 

240  +  'EACH  LINE  OF  THE  DATA  FILE  MUST  BEGIN  * 

250  4-  'WITH  A  LINE  NUMBER' /'FOLLOWED  BY  A  COMMA. '//' EACH  ENTRY  ' 

260  +  'EXCEPT  THE  LAST  ON  ANY  LINE  SHOULD  BE  FOLLOWED  BY  A  COMMA.’// 
270  4-  'FOLLOW  THE  ABOVE  DATA  WITH  LINES  CONTAINING’/'( IN  ADDITION' 

28C  +  '  TO  THE  LINE  NUMBER):'/ 

290  +  '  A)  THE  NUMBER  OF  STORAGE  TEMPERATURES.  (ZERO  IF  NONE.)'/ 

300  +  '  *IGNORE  STEPS  B,C,D,AND  E  IF  THE  NUMBER  ABOVE  IS  0.**'/ 

310  +  ’  B)  THE  PARTICULAR  STORAGE  TEMPERATURES  (DEGREES  KELVIN)'/ 

320  +  *  C)  THE  NUMBER  OF  VALUES  FOLLOWING  (ZERO  IF  NONE.)'/ 

330  +  ’  D)  EITHER  Cl/CO  OR  TIME,  WHICHEVER  IS  TO  BE  SOLVED  FOR.’/ 

340  +  '  E)  PARAMETER  VALUES  FOR  TIME  (YEARS),  OR  Cl/CO  (PERCENT)’ 

350  4-  /’  -DEPENDING  UPON  WHAT  WORD  WAS  TYPED  FOR  PART  D.’//) 

360  STOP 

370  40  WRITE(9,50) 

380  50  FORMAT(/'NAME  OF  INPUT  DATA  FILE’) 

390  READ(9,20)NAM 
400  CALL  DEFINEd.NAM) 

410  WRITE(9,70) 

420  70  FORMAT(/ ' NAME  OF  OUTPUT  DATA  FILE') 

430  READ(9,20) NAM 
440  CALL  DEFI N£(2,NAM) 

450  LN0=0 

460  READ( 1 ,80)LNS,M 
470  80  F0RMAT(2I6) 

480  S3  :0  .  DO 
490  WRITE(9,90) 

500  90  FORMA T(// 'INPUT  DATA’/’  PHI  TEMP '9X ' 1/T ’ 1 1 X 'PHI/TSQUARE ' ) 
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510  L:0 

520  DO  130  l-l  , M 

530  READC 1 , 100) LNS, VEE ,U 

5  AO  100  FORMA  T  C  !  5  ,3010.0) 

550  DO  130  J= 1 ,7 

560  IFCUCJ))130,130,110 
570  1 10  L:L+1 
580  EXCL):1 ./UC J) 

5S0  WYE<L)=VEE/CU(J)**2) 

600  YOU : U  C  J ) 

6  1  0  DVEE=VEE 
620  X  C  L  »  l  )  -  1  • DO 
630  XC  L, 2) : 1 .DO /YOU 

6 AC  WRITE (9, 120) VEE ,UC J> ,EXCL),WYtCL) 

650  120  FORMA T( F 7 . 2  ,  Fg  .2 , 2E  I 4.4) 

660  YCL):DLOGCDUEE/C  YOU* *2 > ) /DLOGC 1  0  .  uO  ) 
670  WYECL):YCL) 

680  S3:S3+YCL.)**2 

690  130  CONTINUE 

700  INO-LNO+l 

710  WH I  I EC 2 , 1 40 ) LNO , L 

720  140  FORMA  I ( 14, 1 H,  ,14) 

730  DO  150  l:l,L 
74C  LNO-LNO+l 

n  I  yw  i  T f  9  KflM  m  t»i  i  \  yv w  r  i  \ 

760  160  F0RMATU4^2(  i’h^EIa!^))  ’  ** 

770  N:2 

7«C)  D0220  I:  I  ,N 
790  DO  200  Jr  1 , N 
800  AC  I ,J) =0.00 
810  DO  180  K= 1 ,L 

820  180  ACI,J):A(I,  D+XCK,I)*XCK,J> 

830  IF  CI.EQ. I ) GOTO  200 
340  LNOrLNO+l 
850  MJ=A;I  ,J> 

860  WHITE  C2, 190)LN0,A1 J 
870  190  FORMA TCI4,lH,,Ei4.6) 

830  200  CONTINUE 

890  9CI)r0.D0 

900  DO  210  K -  1 , L 

910  210  B(I):BCI)+X<X,I)*YCK> 

920  LN0=IN0+1 
930  BI=BCI> 

940  220  WRITE  ( 2 , 1 90 ) LNO ,B I 
950  LNO: LNO+ 1 
960  ESS: S3 

970  WRITE  (2,190) LNO, ESS 
930  DO  230  K: 1 ,N 
990  PrACK.K) 

1000  DO  230  J : K , N 

1010  230  ACK, J):ACK,J)  'P 

1020  8(K):DCK)/P 

1030  DO  250  J : 1 ,N 

1040  IF  CI.EQ. K)  GO  TO  250 

1050  F:ACI,K) 
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1060  DO  240  J=K,N 

1070  240  A(I,J)=A(I,J)-F*A(K,J) 

1080  B(I)=B(I)-F*B<K) 

1090  250  CONTINUE 
1100  WRITE<9,260) 

1110  260  FORMAT <  ///’LEAST  SQUARES  SOLUTION  COEFFICIENTS  ARE:’/) 
1120  DO  280  1  =  1, N 
1130  IMl=I- 1 
1140  BI-B(I) 

1150  WRITE(9,270)IM1,BI 

1160  270  FORMATC  ’  AC, 13,*)  =  '.E15.7) 

1170  LN0=LN0+1 

1180  280  WRITE  <2, 190)LN0,BI 
1190  WRITE(9,290) 

1200  290  FORMAT (//'POLYNOMIAL  CURVE  FIT  OF  DEGREE  l  // 

1210  +  3X 'X-ACTUAL' 7X'Y* ACT UAL' 7X 'Y* CALC  *9X  *DIFF*9X 'PCT-DIFF  /) 

1220  S1=0.D0 

1230  S2=0.D0 

1240  DO  330  1=1  ,L 

1250  C<I)=B(N)*X(I,2)+B(N-1> 

1260  D=Y<I)-C(I) 

1270  SI =S1+D 
1280  S2=S2+D*D 
1290  CEE=C( I ) 

1300  DEE=D 

1310  IF<  C(I).EQ.O.DO)  GO  TO  310 
1320  0bT=100.D0*D/C(I> 

1330  WRITE(9f300)EX(I),WYE(I),CEE,DEE,0UT 
1340  300  FORMATCE 13. 6, 3E1 5.6 ,F1 1 • 1 ) 

1350  G0T0330 

1360  310  WRITEC9 ,320)ZX( I ) ,WYE(I) ,CEE, DEE, RATIO 
1370  320  F0RMAT(E13.6,3E15.6,3X,4A2) 

1380  330  CONTINUE 
1390  EL=L 
1400  EST=S 1/EL 

1420  340TF0RMAT(//'MEAN  ERROR  OF  ESTIMATE  FOR  Y  =  *,E14.6) 

1430  IF(  L-2.NE.0)G0T0  360 

1450  350 ^FORMAT 'NUMBER  OF  POINTS  EQUALS  DEGR|^115JJJUmc°/E  / 
1460  +  'CURVE  FIT  IS  SOLTUION  OF  SIMULTANEOUS  EQUATIONS  / 

1470  +  ’STANDARD  ERROR  DOES  NOT  APPLY') 

1480  GO  TO  380 

1490  360  DF=L-2 

1500  STERR=DSQRT(S2/DF) 

1520  370TF0RMATCUnIiASED  STANDARD  ERROR  OF' /'ESTIMATE  FOR  Ys 
1530  +,16X,E14.6) 

1540  LN0=LN0+! 

1550  ESS=S2 

1560  WRITEC2, 190)LN0,ESS 
1570  380  R= 1 .987D0 
1580  E0=-DL0G( 10.D0 )*R*B(2)*l .D-3 
1590  E=-DL0G< 10.D0)*B(2) 

1600  WRITE(9,390) 
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1610  390  F0RMAT(/’CALCULATED’/3X, 'l/T  ’ 7X ’PHI/TSQUARE ’ ) 

1620  DO  400  1=1  ,L 

1630  F=DEXP(C(I )*DLOG( 10.D0) ) 

1640  IF(I.EQ.L)EFF=F 
1650  FF=F 

1660  400  WRITE(9,300)£X(I) ,FF 
1670  D=DEXP(E*X(L,2))*EFF*E/6G.D0 
1680  EEO  =  EO 
1690  DE£=D 

1700  VRITE(9,410)EE0,DEE 

1710  410  FORMAT(//’E  (KCAL/MOLE) 'E14.6/’A  (SEC-1)  ’,E13.6) 

1720  READ(1,80)LNS,M0 

1730  IF(MO.EQ.O)GOTO  580 

1740  READC 1,100)LNS,(TCI), 1=1 ,M0) 

1750  DO  420  1=1, MO 
1760  TEE=T(I ) 

1770  AK(I)=D*DEXP(-E/TEE) 

1780  420  AKI(I)=AK(I) 

1790  READC 1 ,80)LNS,M1 
1800  I F ( M 1 .  EQ  •  0 )  GOT 05  40 
1810  READC l, 430)LNS, VARY 
1820  430  FORMATC 15 ,2A2) 

1830  IFC VARY. NE. ’TIME’ )WRITE(9, 440) 

1840  440  FORMATC /’DEGREES  K’8X’PER  CENT,  Cl  OF  CO,  (TIME  IN  YEARS)’) 

1850  IFC VARY. EQ. ’TIME’ ) WRITE C 9, 460) 

1860  460  FORMATC/ ’DEGREES  K*8X*TIME  IN  MINUTES,  CC1/C0  IN  PERCENT)’) 
1870  470  FORMATC ’KELVIN  (SEC  -l)’,6F8.1/> 

1880  READ  C 1 , 100) LNS, (Z(J),J=1,M1) 

1390  DO  490  J=l,Ml 
1900  ZEE=Z(J) 

1910  IFC VARY. NE. ’TIME* )  Q( J)=ZEE*365.DO*24.DO*3SOO.DO 
1920  IFC VARY. EQ. ’TIME’ )  QC J)=DL0GC 100.D0/ZEE)/60.D0 
1930  490  CONTINUE 
1940  WRIIE(9,470)(ZCJ),J=1,M1) 

1950  DO  520  1=1, MO 
1960  DO  510  J= 1 ,M1 
1970  IF(VARY.EQ.’TIME')G0T05O0 
1980  F=Q( J)*AK(I) 

1990  IF(F.GE.9.91)DEG(J)=0. 

2000  IF(F.LT.9.91D0)DEGCJ)  = 100 „DO*DEXP( -F) 

2010  GO  TO  510 

2020  500  DEGC J) =QC  J)/AK(  I ) 

2030  510  CONTINUE 

2040  520  WRITE(9,530)T(I),AKI(I),(DEG(J),J=1,M1) 

2050  530  F0RMATCF7.1 ,E10.2,2X,6F3.2) 

2060  GO  TO  580 
2070  540  WRITEC9 ,550) 

2080  550  FORMATC ’DEGREES  X’/’KELVIN  CSEC-1)’) 

2090  WRITE(9,560)(T(I),AKI(I),I=1,M0) 

2100  560  FORMATC 4CF7. 1 ,E10.2) ) 

21 10  580  VRITEC9,590) 

2120  590  FOKM,.T(///’S-  OUTPUT  FILE’) 

2130  STOP 
2140  END 
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DIMENSION  NANC6) 

COMMON  X<30> ,YC30) ,N,EN,SXC2) , SY<2) ,SXYC2> ,A,BC2) ,XBAR,YBAR 
DOUBLE  PRECISION  X,Y,EN,SX,SY,SXY,A,B,XBAR,YBAR 
WRITEC9, 10) 

INPUT  DATA  FILE* ) 


10  FORMATC 'NAME  OF 
R£ADC9,20) NAM 
20  FORMAT (6A2) 

CALL  DEFINEC 1 ,NAM) 
WRITEC9,30) 

30  FORMATC/.  INPUT 
100  READC 1 ,40)LN0,N 
110  40  F0RMATC2I6) 


DATA  VAX  *X  '  ,  13X  *Y  * ) 


170  DO  SO  1=1, N 

180  READC1,50)LN0,X(I),YCI) 

190  50  FORMATC 16 ,3D20 .0) 

191  EX=XCI) 

192  EY=YCI) 

193  SO  WRITEC9,70)EX,EY 
200  READC 1,50)LN0,SYC1> 

205  READC 1 ,50) LNO,SXC 1 ) 

210  READC 1,50)LN0,SXC2) 

215  READC 1 ,50)LN0,SXYC 1 ) 

220  READC 1 ,50)LN0,SYC2) 

225  READC  1 ,50)LN0, BCD 
230  READC 1 ,50) LNO,B C2> 

235  READC 1 ,5O)LN0,SXYC2) 

260  70  F0RMATC2E14.6) 

270  EN=N 


280  YBAR=SYCl)/EN 

290  A=YBAR 

300  XBARrSXC 1 )/EN 

310  SXYC2)=SXYC2)/CEN-2.D0) 

312  SA=A 
314  SB=BC2) 

316  SXBAR=XBAR 

320  WRITEC9,80)SA,SB,SXBAR 

330  80  FORMATC/ *Y=A+BCX“XBAR) */*Ai*,El4.6/*B|* ,E14.6/*XBAR*,E14.6) 

341  IFRAT=1 

342  90  FORMATC 16 ,2G10 .0) 

350  DO  95  1=1,2 

352  READC 1 ,90)LN0, CL, F2 

360  95  CALL  INDEPCI, CL, F2,IFRAT) 

361  DO  110  1=1,2 

362  READC 1 , 100)LNO,CL,T ,CHISQ1 ,CHISQ2 
364  100  FORMATC 16, G10. 0,3010.0) 

370  110  CALL  CONLIMC I ,CL,T,CHISQ1 ,CHISQ2) 

380  STOP 
390  END 

1310  SUBROUTINE  INDEPCI ,CL,F ,IFRAT) 

1320  COMMON  XC30) ,YC30) ,N,EN,SXC2) ,SYC2> ,SXYC2> ,A,BC2> ,XBAR,YBAR 
1325  DOUBLE  PRECISION  X,Y,EN,SX,SY,SXY,A,B,XBAR,YBAR 


-22- 


PROG 1  -  PAGE  2 


1530 

1335 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1403 

1406 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1565 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1650 

1655 

1656 

1657 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 


DIMENSION  SQS(3),IDF<3) ,DF<3) ,SQM(3) 

G0T0C5 ,35) , I 

5  SQS< 1 )=<SXY< 1 )-SX<l )*SY< 1 )/EN)**2/(SX<2)-SXC 1 >**2/EN) 

IDF< 1 )= 1 
DF(1)=IDF(1) 

SQM(1)=SQS(1)/DF(1) 

SQS(3)=SY(2)-SY(1)**2/EN 
IDFC3)=N-1 
SQS(2)=0.D0 
DO  7  IH0PE=1,N 

7  SQS(2)=SQS(2)+(Y(IH0PE)-(A+B(2)*(X(IH0PE)-XBAR)) )**2 
IDF(2)=IDF(3)-IDF(l) 

DF(2)=IDF(2) 

SQM(2)=SQS(2)/DF(2) 

FRATI0=SQM(1)/SQM<2) 

WRITE(9,10)SQS(1),IDF( 1 ) ,SQKC1)  ,FRATIO 

10  FORMA T(/28HTEST  FOR  SLOPE  EQUAL  TO  ZER0///38X, 10HDEGREES  OF/ 
+23X49HSUM  OF  SQUARES  FREEDOM  MEAN  SQUARE  F  RATIO/ 
+17HDUE  TO  REGRESSION  5X  £14.6, 17 ,5X,E13.3,3X,F7.2) 

WRITEC9 ,20)SQS(2) ,IDF(2),SQM(2) 

20  F0RMATC8KRESIDUAL1 4XE 14. 6,17, 5X, £13. 3) 
WRITE(9,30)SQS<3),IDF<3) 

30  F0RMAT(/5HT0TAL17XE14.6,I7) 

35  WRIT£(9,40)CL,IDF(1),IDFC2),F 

40  F0RMAT(/2HATF5.1,12H  Z  LEVEL  AND, 13 , 1 H,  ,  13 , 1 1 H  DEGREES  OF 
+12H  FREEDOM,  F=,F8.3) 

IF(FRATIO.GT.F)G0T060 

IFRAT=IFRAT+1 

WRITE(9,50) 

50  F0RMATC40HF  RATIO  LESS  THAN  F,  CONSEQUENTLY  CANNOT) 

60  WRITE(9,70) 

70  FORMA T( 47 HACCEPT  THE  HYPOTHESIS  THAT  Y  IS  DEPENDENT  ON  X.) 

RETURN 

END 


SUBROUTINE  CONLIM(J,CL,T,CHISQi ,CHISQ2> 

COMMON  X<30; ,Y(30),N,EN,SX(2) ,SY(2) ,SXY(2) ,A,B(2) ,XBAR,YBAR 
DOUBLE  PRECISION  X,Y,EN,SX,SY,SXY,A,B,XBAR,YBAR,VARC3) , 
+SDEV(3) , T ,CHISQ 1 ,CHISQ2 , ACOEFF ,BC0EFF,XC0EFF,SX2,YC 
G0T0C5 , 1 55 ) , J 

5  VAR(l)  =  <SX<2'>-SX<l)**2/EN)/(EN-l.D0) 

VAR<2)=(SY(2)-SY< l )**2/EN)/(EN-l .DO) 

VAR(3)=(EN-1 .DO ) *( VAR (2) -B(2)*B(2)*VAR ( 1 ) )/(EN-2.D0) 

DO  10  1=1 ,3 

10  SDEV( I ) =DSQRT( VAR( I ) ) 

VRITE(9,120)N 


- • - ' ^ I  • / nwm»  fwvwvv  i  /  | v  n  n  V  J  J 

WRITE (9, 1 40)YBAR,SDEV(r i ,VARC2) 

WRITEC9, 150)SDEV(3) ,VAK(3) 

120  FORMATC/llHSAMPLE  SIZE, 14// 1 2HESTIMATES  0F/23X4KMEAN , 
+5HSIGMA ,9X, 13HSIGMA  SQUARED)  * 
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1770  130  F0RMAT(20HINDEPENDENT  VARIABLE, 3D  1 4.6) 

1730  140  F0RMAT< 18HDEPENDENT  VARIABLE, 2X, 3D  1  4.6) 

1790  150  F0RMAT(3HY.X,31X,2D14.6) 

1800  155  AC0£FF=T*SDEV(3)/DSQRT(EN> 

1810  AL0V=YBAR-AC0EFF 
1820  AHI3H=YBAR+AC0EFF 
1325  YB=Y8AR 

1330  V1RITE(9 , 1 60  )  CL,  ALOW,  YB,  A  HIGH 

1340  160  F0RMAT(/F4.l,23H  PER  CENT  CONFIDENCE  LI MI TS://9HPARAMETER, 
1850  +13X,  1  IHLOWER  LIMIT  ,3X,9KF.STIMATED5Xl  1 HUPPER  LIMIT,/ 

1860  +16HI NTERCEPT  (YBAR) ,3X,3E14.6) 

1370  BC0EFF=T*SDEV<3)/CSDEV<1)*DSQRT(EN-1.D0)) 

1330  BL0W=B(2)-BC0EFF 
1890  BHIGH=BC2)+BC0EFF 
1892  BEE:B(2) 

1900  WRITE(9,170)BL0W,B£E,8HIGH 
1910  1 70  FORMA T(5HSL0PE, 1 4X ,3E 1 4.6) 

1920  SYX2L0= VARC3) /CHISQ 1 
1930  SYX2HI=VAR(3)/CHISQ2 
1932  VARYX=VAR(3) 

1940  WRITE(9, 180)SYX2L0,VARYX,SYX2HI 

1950  180  FORMA T( 1 9HSI GMA  SQUARED  (Y.X)3E14.6) 

1960  WRITE(9, 1  90)CL 

1970  190F0RMAK/F5.1  ,41H  7.  CONFIDENCE  LIMITS  FOR  MEAN  OF  Y  VALUES 
1930  +24H  AT  A  PARTICULAR  X  VALUE//4X, 1HX, 13X,6HY  CALC,SX,7HY  UPPER 
1990  +7X,7HY  LOWER) 

2000  AC0EFF=T*SDEVC3) 

2100  SX2=SX(2)-SX(1)**2/EN 
2110  DO  200  1=1, N 

2120  XCOEFF=X(I)-XBAR 

2121  YC=YBAR+BC2)*XC0EFF 
C130  YC 1 =DEXP (YC*DLOG( 10 .DO ) ) 

2140  XCOEFF=DSQRT ( 1 .DO/EN  +  XC0EFF**2/SX2) 

2150  YC2=DEXP( < YC+ACOEFF*XCOEFF)*DLOGC 10. DO) ) 

2160  YC3=DEXP( ( YC-AC0EFF*XC0EFF) *DL0G( 10. DO)) 

2170  EX=X(I) 

2180  200  WRI TE(9 ,210 )EX,YC 1 , YC2 ,YC3 
2190  210  F0RMATC4E14.6) 

2200  RETURN 
2210  END 
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10  DIMENSION  NAMC6) ,NC3) , IFHATC2) , ITRATC2,2) 

?0  DOUBLE  PRECISION  EN(  3  )  ,X<  2,  50  ) ,  Y(  2 , 50 )  ,SX  (3 ,2 )  ,SY  ( 3 , 2 ) 

30  f,SXYC3,2> ,XBAfiC3),YBAR(3),\/ARX<3) , VARY <3 )  ,B<3 f 2) , VA RXYC3 ) 

AO  DO  10  1:1 ,3 
50  N C I ) : 0 
(SO  DO  10  J:  1  ,2 
70  SX<I,J):C. 

30  SY< I, J)=0. 

90  10  SXYCI,J):0. 

100  I S: 1 

110  20  WR l TE C  9 ,30)  I S 

120  30  FORMA T( /23H NAME  OF  INPU1  DATA  FILE, 12) 

130  REA  DC  9 , AO ) NAM 
1  AC  AO  FORMA T ( 6A2 ) 

150  CALL  D£F I  NEC  IS , NAM) 

160  JR  I TE ( 9 , 50 ) 

170  50  F0RMATC//1 1HINPUT  DATA ! /3X , 1 1 M 1 N DEPENDENT , 5 X , 9H DEPENDEN T ) 
180  READC I S , 60 ) LNO , N C IS ) 

190  60  F0RMATC2I6) 

2CC  NS^NCIS) 

210  DO  70  i  = 1  ,NS 

220  REA  DC  I S , 80 ) LNO,X< IS , I ) , YC IS ,  I ) 

230  80  FORMA T ( 16 ,3020.° ) 

2A0  EXrXC IS, I) 

250  £Y  =  Y(  13, 1 ) 

260  70  WR I TE C9 ,90 ) EX , EY 
270  R£AD(IS,30)LN0,SYCIS, 1 ) 

280  READCIS,80)LN0,SXCIS,1  ) 

290  READC I S , 30 ) LNO, SXC IS ,2 ) 

300  READC IS, 80 ) LNO, SXY C IS , 1 ) 

310  READ(IS,30)LN0,SY(IS,2) 

320  REA DCIS, 80)LN0,B(IS,1 ) 

330  READCIS,80)LN0,B(I3,2) 

3 AC  READC IS,80)LN0,SXYC IS ,2) 

350  90  FORMAT C2E 1  A, 6) 

360  I S  =  IS-*- 1 

370  IF( IS-2)20,20,IC0 

380  100  DO  1  AO  1  =  1 , IS 

390  GOTOC 1 10,1 10,130)  , I 

ACC  1  1C  DO  120  J=1  ,2 

A10  SXCIS,I):SXCIS, I)+SXCJ,I) 

A2C  120  SY<IS,I)=SYC1S,I)+SYCJ,I) 

A30  N(IS)=NCIS)+NC I ) 

4 AO  SXYCIS.l )=SXYC IS, I )+SXYCi ,1 ) 

450  130  ENC I >=NC I ) 

A60  XBARCI )=SX( I ,  I  )/ENCI) 

470  YBARC I ) =SYC I , 1 )/£N(J) 

480  WARXCI)=SXCI,2)-SX(1,I )**2/ENCI) 

490  VARYC I )  =  SY C I ,2)-SYC I , 1 )*+2/EN( I ) 

500  SXYC  I,2):SXYC I , 1 ) -SXC I , 1 )*SYf 1,1 )/ENCl ) 


PR0G2  -  PAGE  2 


510  B!IS,2)=SXY!I,2>/VARX!I) 

520  1 40 VARXY !  I)=!VARY!I) -B(  1 S  ,  2  )**2*  VAHX  C 1 ) )  /  ( LYl  1 )  -2.  DC 
530  EFF:VARXY! 1 )/VAHXY!2) 

540  IDF) =  N! l )-2 
550  IDF2=N!2)-2 
560  WR1TE!9, 150HDF1  ,IDF2,EFF 

57C  15CFORMAT!/I4,lH,,l4, *  DEGREES  OF  FREEDOM ’ / 1 1 0 . 4  ,  *  F  RATIO’//) 

571  DO  ISC  1:1,2 

572  READ  !2,2OO)LN0,CL,F 

573  I FRATC 13=1 

574  WRITE  (9  ,  l  60 )C  L »r 

575  160F0RMAKF7. 1  , 19H  7.  CONFIDENCE  LEVEL/ 

576  +.F7.3.24H  VALUE  OF  F  DISTRIBUTION) 

577  IF(EFF.L£.F)GOT0130 
573  I FRaT< I ) :2 

57?  WR I T£!9 , 170) 

530  170  FORMAT! 1 9riS I GN I F I CA NT ,  CANNOT  ) 

581  ISO  WR 1 TE ( 9 , 1 90 ) 

582  190  FORMAT ( 4 4H ACCEPT  THE  HYPOTHESIS  THAT  THE  VARIANCES  ARE 
533  +6 H  EQUmL/) 

584  I F( ( 1  FRA T ( 1 ) . Ed . 2 ) . AND . (I FR AT ( 2 ) . EQ . 2 ) ) G0T035O 
535  31 «u2:3Cl ,2)-SC2,2> 

590  r>F:N(  1  )+N<2)-4 

591  NDFrDF 

600  SYXP2=(!EN! 1 )-2.D0)* VARXY! l )+! EN ! 2) -2. DO ) * VARXY! 2) )/DF 
610  VARB71B:SYXF2*(  I .  DO/VARX(  1  )  + 1 . DO /VARX  ( 2 ) ) 

620  TEE:Bli’®2/S3RT(VARBMB) 

621  DO  290  1=1 ,2 

630  READ(2,200)LN0,CL,T 

635  2CC  FORMAT!  16,3010.0) 

636  I TRA T< 1,15  =  1 

637  IF( I.EQ.2) G0T022O 

638  WRITE(9,21C)NDF,TEE 

639  210  FORMA T C / 17 , 1 SH  DEGREES  OF  FREEDOM/ 

640  +F7.3.22H  VALUE  OF  T  CALCULATED/) 

641  220  IF<IFRAT(I).E(i.2)GOT0290 

642  IF! ABS!TEE) .GT.T)  ITRAT! 1 ,19-2 
720  VRITE(9,230)CL,T 

730  230  FORMA T C/F7 . 1  , 1  9H  7.  CONFIDENCE  LEVEL 

740  +/F7.3,33rl  ABSOLUTE  VALUE  OF  T  DISTRIBUTION.' 

741  J  F( I TRAT! 1 , 1 ) ,£Q . 2) GOT 0260 

742  WRITE! 9,250 ) 

743  250  FORMAT! 40HNON-SIGWI FiCANT  ! THE  LINES  ARE  PARALLEL)//) 

744  GOT0290 

745  260  WR I TE ( 9 , 270  > 

746  270  FORMAT!  /40HSIGNIFICANT  (THE  LINES  ARc.  NOT  PARALLEL)//) 

74/  290  CONTINUE 

74  3  IF! ! I  THAT! I  ,1 ).EQ,2) .AND.! I  THAT! 1 ,2) .EQ.2) ) G0T03  5C 
780  XVAriRL:  VARX!  1  )  +  VARX!2) 

790  8(3,2)  =  ! SXY ! 1 , 2 )+SXY! 2 ,2) ) /XVARPL 

80 C  SXY T VO: VARY! 1 )+ VARY ! 2 ) -Bi3 ,2 ) * ! SXY ! 1 ,2)+SXY!2,2) ) 
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810  DP: N( 1 )+N( 2 ) ~3 

820  NDF:DF 

330  SYXP2=SXYTW0/DF 

840  XBARMI=X8AR( 1)-XBAR(2) 

850  TEE=(YBAR(  1>-YBAR(2)-B(3,2)*X6ARMI)/S(IRT(SYXP2*(  1  .DO/ENC  1  >  + 
360  +1.D0/EN(2)  +  XBARfll **2/X VARPL) ) 

86 l  DO  340  1=1,2 

370  READ<2,200)LN0,CL,T 

871  I TRA  T(2 , 1 )  =  I 

372  IF< I . EQ .2 ) GOT 0295 

373  WRITE(9,210)NDF,TEE 

874  295  I F( ( I FRATC l).EQ.2).0R.( ITRAT ( 1 , 1 ) . £3.2) ) G0T03 40 
37  5  1F(A9S(TEE) ;0T: T)1TRAT(2,I )  :2 

876  4RITE(9,23C)CL,T 

877  ]F( ITRAT(2,I).EQ.2)G0T0320 
900  WRI TE ( 9 , 3 1 0  ) 

910  3!0F0RMAT(/4?H NGN “SIGNIFICANT  (THE  LINES  CAN  B £  REGARDED  AS  A 
920  +24H  SINGLE  COINCIDENT  LINE)) 

930  G0T0340 

940  320  WR ITE<  9  ,330  ) 

950  330FORMAT( 42HSIGN1FICANT  (THE  LINES  ARE  NOT  COINCIDENT)) 

970  340  CONTINUE 
980  350  STOP 
990  END 
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Preceding  page  blank 


M-26  CONTROL 


INPUT  DATA 


PHI 

TEMP 

1/T 

PHI/TSQUARE 

80,00 

497.57 

0.2010E-02 

0.3231 E-03 

60.00 

498.67 

0 .2005E-C2 

0.3217E-03 

80  o  00 

497.67 

0.2009E-02 

0.3230E-03 

80.00 

497,48 

0 . 20  i  OE  -02 

0.3233E-03 

80.00 

490.72 

0.2038E-02 

O.3322E-03 

40.00 

486.35 

0 , 2056E-02 

0. 1691E-03 

40.00 

486.55 

0 . 2055E-G2 

0, 1690 E-03 

40.00 

490,21 

0.2C40E-02 

0. 1665E-03 

40. OC 

483.00 

0.2070E-02 

0. 1715E-0S 

20.00 

463.27 

0  4  2069E  -02 

0 .8563E-04 

20.00 

481.70 

0.2076E-02 

0.3619E-04 

20.00 

482.58 

0 . 2072E -02 

0. 8588E-04 

20.00 

477.36 

0.2093E-02 

0 . 3V58E-04 

10. GO 

A  T  C  AT 
*1  1  «/  •  V  • 

C. 21052-02 

C.4431E-04 

10.00 

476.25 

0.2100E-02 

0 . 4409E-04 

10.00 

478.11 

0.2092E-02 

0.4375E-04 

1C. 00 

476.34 

0.2099E-02 

0.4407E-04 

5.00 

471.39 

0.2121E-02 

0 . 2250E -0  4 

5. GO 

471,29 

0.2122E-02 

0.2251E-04 

5.00 

474.91 

0.2106E-C2 

0 .221 7E-04 

5,00 

469.24 

G.2131E-C2 

0.2271E-04 

2.50 

466.42 

0.21 44E-02 

0 . 1 1 49E-04 

2.50 

465.25 

0.21 49E-02 

0. 1 l 55E-04 

2.50 

465.44 

0.2149E-02 

0.1 I54E-04 

2.50 

464.67 

0,21 52E-02 

0. 1 1582-04 

E  (KCAL/MOLE)  0.492896E  02 
A  (SEC-1 )  0.748259E  21 

DEGREES  K  PERCENT,  UNREACTED,  (TIME  IN  YEARS) 


KELVIN 

(SEC  -l  ) 

1  .0 

2.0 

5.0 

10.0 

20.0 

298.0 

0.53E-15 

100.00 

100.00 

too.oc 

100.00 

100. CO 

323.0 

0.33E-  12 

IOC. 00 

100.00 

99.99 

99.99 

99.98 

3  43,0 

0.33E-10 

99.74 

99.48 

98.71 

97.43 

94.93 

373.0 

0.98E-03 

73.39 

53.36 

21 .29 

4.53 

0.21 

-30- 


M-26  /  MK  12  GUN  PLUG  FOAM 


INPUT  DATA 


PHI 

TEMP 

1/T 

PHI/TSQUARE 

80.00 

494.99 

0.2020E-02 

0.3265E-03 

80.00 

495. 6S 

0.2017E-02 

0.3256E-03 

40.00 

484.87 

0.2062E-02 

0.1 701E-03 

40.00 

433.49 

0.2068E-02 

0.171  IE-03 

40.00 

480.14 

0.2083E-02 

0. 1735E-03 

40.00 

485.07 

0.2062E-02 

0. 1700E-03 

20.00 

483 .86 

0.2067E-02 

0.3543E-04 

20.00 

433.67 

0.206SE-02 

0.8549E-04 

10.00 

476.93 

0 . 2097  £-02 

C.4396E-04 

10.00 

475.66 

0.2102E-02 

0.4420E-04 

5.00 

471 .39 

0.2121 E-02 

0, 2250E-0  4 

5 .00 

4vS . 63 

n  oi 

v  •  &  i  Oc* 

C . 22S7E-C  4 

2.50 

468.17 

0.2136E-02 

0.1 1 41E-04 

2.50 

464.47 

0.21 53E-02 

0. 1 159E-04 

E  (KCAL/MOLE) 

0. 

533624E  02 

A  (SEC- 

1  ) 

0. 

l  12227E  24 

DEGREES 

K 

PERCENT, 

UNREACTED,  (TIME  IN  YEARS) 

KELVIN 

(SEC 

-1 ) 

l  .0 

2.0 

5.0 

10,0 

20.0 

298.0 

O.OOE 

00 

100.00 

100.00 

100.00 

100.00 

100.00 

323.0 

0.40E 

-13 

190.00 

100.00 

100,00 

100.00 

100.00 

343.0 

0. 17E 

-10 

99.95 

99.90 

99.74 

99.48 

98.96 

373.0 

0.31E 

-03 

90.75 

82.36 

61.55 

37.89 

14.36 

-32- 
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M-26  /  MK  12  GUN  PLUG  FOAM, 
HIMIDITY-AGED 


;nput  data 

phi 

TEMP 

1/T 

PHl/ISQUARE 

40.00 

431.08 

0.2079E-02 

0. 1 728E-03 

40.00 

481.57 

0.2077E-02 

0. 1725E-03 

20.00 

434,44 

0.20S4E-02 

0.8522E-04 

20.00 

482.95 

0.207IE-02 

0.8575E-O4 

10.00 

477.95 

0.2092E-02 

0.4573E-C4 

10.00 

474.50 

0.2107E-02 

0.444IE-04 

5.00 

467.67 

0.2I38E-02 

0 . 2286E-0  4 

5.00 

470, SO 

0 . 2 1  25E-02 

0.2258E-04 

2.50 

467.J4 

0.2I40E-02 

0.11 45E-04 

2.50 

465.95 

0.2155E-02 

0.1 162E-04 

1  .25 

459.52 

0.2177E-02 

0.5925E-05 

1.25 

458.55 

A  n  «  r\  r* 

\j  •  1  os.c  -u  c 

0 , 5950E-05 

E  (KCAL/MOLE)  0.5458S9E  02 
A  (SEC-l)  0.265494E  24 

DEGREES  K  PERCENT,  UNREACTED,  (TIME  IN  YEARS) 


KELVIN 

(SEC  -1 ) 

1  .0 

2.0 

5.0 

10.0 

20.0 

293.0 

C.OOE  00 

100.00 

100.00 

100.00 

100.00 

100,00 

323.0 

C.30E-13 

100.00 

100.00 

100.00 

100.00 

100.00 

348.0 

0. 14E-10 

99.96 

99.91 

99.78 

99.57 

99.14 

373.0 

0.27E-03 

91  .76 

84.19 

65.04 

42.30 

17.90 

M-26  /  LOCTITE  A 


INPUT  DATA 


PHI 

TEMP 

l/T 

PHI/TSQUARE 

80.00 

493.50 

0 . 2026E-C2 

0.3285E-03 

40.00 

488.53 

0.2047E-02 

0. 1 676E-03 

20.00 

479.73 

0.20S5E-02 

0 . 8690E-0  4 

10.00 

473.12 

0.211 4E-02 

0. 4467E-04 

5.00 

469.53 

0.21 30E-02 

0.2268E-04 

2.50 

463.98 

0,21 55E-02 

0. 1  16  IE-04 

E  (KCAL/MOLE)  0.497688E  02 
A  (SEC-1)  0.150313E  22 


DEGREES  K  PERCENT,  UNREACTED,  (TIME  IN  YEARS) 


KELVIN 

(SEC  -1  ) 

l  .0 

2.0 

5.0 

10. 0 

20.0 

298.0 

0.47E-15 

100.00 

IOC, 00 

100.00 

100.00 

100.00 

323.0 

U.32E-1 2 

100, cc 

100.00 

100.00 

99.99 

99.98 

348.0 

0.S3E-1O 

99.74 

99.48 

98.70 

97.42 

94.90 

373.0 

0.  10E-G7 

72.21 

52.  1  5 

19.64 

3.36 

0.15 

-:*G- 


A _ 

.0019 


Act  tv  At  ion 

jttEQUElijCY  FACTOR  *  .15  X 
pTATISlriCALL|r  COMPATIBLE 


4918  KCAL  per  HOLE 


1022  per  iECOND 


J . 


.0020 


1/T 


.0021 


Best  Available  Copy 
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M-26  /  LOCTITE  N 


INPUT  DATA 


Pril 

TEMP 

1/T 

PH1/TSQUARE 

80.00 

497.97 

0.2008E-02 

0.3226E-03 

40.00 

483.73 

0.2046E-02 

0. 1 675E-03 

20.00 

479.34 

0.2086E-C2 

0.8704E-04 

20.00 

430. 42 

0.20S2E-02 

0 . 8665E-04 

10.  OC 

476.15 

O.2IOOE-02 

0.441 IE-04 

5.00 

471.98 

0.21 19E-02 

C.2245E-04 

2.50 

464.57 

0.2153E-02 

0. 1 158E-04 

E  (KCAL/NOLt)  0.474187E  02 
A  ( SEC- l )  0 . 1 03894E  21 

DEGREES  K  PERCENT,  UNREACTED,  (TIME  IN  YEARS) 


KELVIN 

(SEC  -1) 

1  .0 

2.0 

5.0 

10.0 

20.0 

293. C 

0. 17E-14 

100,00 

.'00.00 

100.00 

100.00 

100.00 

323.0 

0.85E-12 

100. CO 

99.99 

99.99 

99.97 

99.95 

343.0 

O.17E-0S 

99.46 

98.92 

97.33 

94.73 

89.74 

373.0 

0. 17E-07 

58.50 

34,22 

6.85 

0.47 

0.00 
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i  /TSQUAAE 


n 

M-26  /  EC  1099 


INPUT  DATA 


PHI 

TEMP 

1/T 

PHI /773UARE 

f  ' 

80.00 

488.04 

0.2049E-02 

0 . 3359E-03 

40.00 

484.77 

0 .206CE-02 

0.1 702E-03 

20.00 

477. OS 

0.209SE-02 

0 . 8787E-0  4 

f 

10. CO 

471 .85 

0.21 19E-02 

0. 4492E-C4 

r- 

5.00 

467.00 

0 . 21 4 i £-02 

0  *  2293E-0  4 

* 

2.50 

464.96 

0.21 5  IE -02 

0.1 1 5SE-0  4 

r 

? 


E  (KCAL/MOLE)  0.596293E  02 

A  (SEC-1 )  0.806879E  26  !  I 

DEGREES  K  PERCENT,  UNREACTED,  (TIME  IN  YEARS) 


KELVIN 

(SEC  -1  ) 

l  .0 

2.0 

5.0 

10.0 

20.0 

293.0 

0.00E  00 

100.00 

100.00 

100.00 

100.00 

100.00 

323.0 

0 .COE  00 

100.00 

100.00 

100.00 

100.00 

100.00 

343.0 

0.29E-11 

99.99 

99.98 

99.96 

99.91 

99.32 

373.0 

0.92E-09 

97.13 

94.34 

86.45 

74.74 

55.86 

Li 


r 


-4U- 


. . nHp  i  rt  wm 


M-26  /  UNFILLED  CARTRIDGE  CASE 


INPUT  DATA 


PHI 

TEMP 

1/T 

PHI /TSQUARE 

80.00 

494.99 

0.2020E-02 

0.3265E-03 

40.00 

488.73 

0 , 2C4SE-02 

0 . 1 675E-03 

20.00 

485.64 

0 . 20  5  9  E  -0  2 

0.8480E-04 

10.00 

476.54 

0.209SE-02 

0 .4404E-04 

5.00 

471.88 

0.21 19E-02 

0.2245E-04 

2.50 

464.47 

0.2153E-02 

0. 1 159E-04 

E  ( KCAL/MOLE )  0.496687E  02 

A  (SEC-1)  0.103975E  22 

DEGREES  K  PERCENT,  UN^EACTED,  (TIME  IN  YEARS) 


KELVIN 

(SEC  -1) 

1  .0 

2.0 

5.0 

10,0 

20. C 

298.0 

0.39E-15 

100.00 

100,00 

100.00 

100.00 

100.00 

323.0 

0.26E-12 

100.00 

100.00 

100.00 

99.99 

99.98 

348.0 

0.66E-10 

99.79 

99.58 

98.96 

97.93 

95.90 

373.0 

0.82E-08 

77.28 

59,72 

27.56 

7.60 

0.58 

M-26  /  FILLED  CARTRIDGE  CASE 


INPUT  DATA 


PHI 

TEMP 

l/T 

PHI/TSQUARE 

80.00 

496.98 

0.2012E-02 

C.3239E-03 

40.00 

491.50 

0.2035E-02 

0 . 1 65SE-03 

20.00 

485.15 

0.2061E-02 

0.8497E-04 

20.00 

484.95 

0.2C62E-02 

0.8504E-04 

5.00 

471.7s 

0.2120E-02 

0.2246E-04 

2.50 

465.44 

0 .21 49E-02 

0. 1 154E-04 

E  (KCAL/MOLE)  0.477363E  02 
A  (SEC-1 )  0.1 1624  IE  21 


DEGREES  K  PERCENT,  UNREACTED,  (TIME  IN  YEARS) 


KELVIN 

(SEC  -1  ) 

1.0 

2.0 

5.0 

10.0 

20.0 

298.0 

0.1  1E-U 

100.00 

100.00 

100,00 

100.00 

100.00 

323.0 

0.58E-12 

100.00 

100.00 

99.99 

99.98 

99.96 

348.0 

0. 12E-09 

99.62 

99.24 

98.11 

96.25 

92.64 

373,0 

0. 12E-07 

67.65 

45.76 

14.17 

2.01 

0.04 

-44- 


PYRO  CONTROL 


INPUT  DA' 


TEMP 


PHI /TSQUARE 


DEGREES 

KELVIN 

298.0 


80.00 

503.47 

0. 1986E-02 

0.31  56E-03 

40.00 

496.90 

0,201 2E-02 

0 . 1 620E-03 

40.00 

498.39 

0.2006E-02 

0. 1610E-03 

40.00 

497.45 

0.2010E-02 

0.161 6E-03 

20.00 

483.93 

0.2045E-02 

0.8366E-04 

20.00 

489.43 

0.2043E-02 

0.3349E-04 

20.00 

489.49 

0.2043E-02 

0.8347E-04 

10.00 

482.86 

0.207  IE -02 

0 . 4289E-0  4 

10.00 

482.86 

0.207 1 £-02 

0 . 4289E-0  4 

10.00 

483.32 

0.2069E-02 

0.4281 E-04 

5.00 

476.98 

0.2097E-02 

0.21981  -0  4 

5.00 

477.43 

0.2094E-02 

0.2193E-04 

5.00 

475.80 

0.2102E-02 

0 . 2209E-0  4 

2.50 

468.81 

0.2133E-02 

0.1  137E-04 

2.50 

468 .32 

0.2135E-G2 

0.1  140E-C4 

2.50 

469.93 

0.2128F-P2 

0.1  132E-04 

1.25 

462.71 

0 .2 1  6  1  E  -02 

0.5838E-05 

1.25 

462.71 

0,2161 E-02 

0.5838E-05 

MOLE) 

0.443645E  02 

) 

0.193012E  19 

K 

PERCENT,  UNREACTED,  (TIME  IN  YEARS) 

(SEC 

-1) 

1.0  2.0 

5.0  10.0 

20. 

0.5SE 

-14  100.00  100.00 

100.00  100.00 

100 

0.13E 

-11 

99.99  99.99 

99.97  99.94 

99 

0.26E 

-09 

99.17  98.35 

95.92  92.01 

84 

0. 19E 

-07 

54.13  29.30 

4.65  0,22 

0 

PYRO  /  UNFILLED  CARTRIDGE  CASE 


INPUT  DATA 

PHI  TEMP  |/T  PHI/TSQUARE 

80.00  503.  97  0.1984E-02  0.3150E-03 

AO. 00  495.91  0.201 6E-02  0.1627E-03 

40.00  494.91  0.2021E-02  0.1633E-03 

40.00  495.47  0.2018E-02  0.1629E-03 

20.00  490.43  0.2039E-02  0.8315E-04 

20.00  488.43  0.2047E-02  0.8333E-04 

20.00  439.00  0.20  15E-02  0.8364E-04 

10.00  431.87  0.2C75E-02  0.4307E-04 

10.00  481.87  0.2075E-02  0.4307E-04 

10.00  481.55  0.2077E-02  0.4312E-04 

5.00  475.51  0 ,2 103E-02  0.22I1E-04 

5.00  475.51  0.2103E-02  0.2211E-04 

5.00  475.11  0 .2 105E-02  0.2215E-04 

2.50  463.31  0.2133E-02  0.1137E-04 

2.50  468.32  0.2135E-02  0.1140E-04 

2.50  469.54  G.2130E-G2  u,li34E-04 

1.25  463.20  0 .21 59E-02  C.5326E-05 

1.25  463.20  0,21 59E-02  0.5326E-05 


E  (KCAL/MOLE)  0.457713E  02 
A  (SEC-1)  0.932722E  19 


DEGREES 

K 

PERCtNT, 

UNREACTED,  (TIME  IN  YEARS) 

KELVIN 

(SEC  -l) 

1 .0 

2.0 

5.0 

10.0 

20.0 

298.0 

0.25E-14 

100.00 

100.00 

100.00 

100.00 

100.00 

323.0 

0.99E-12 

.00.00 

99.99 

99.98 

99.97 

99.  S  'i 

348.0 

0.17E-09 

99.48 

98.95 

97.40 

94.88 

90.01 

373.0 

0, 14E-07 

64.12 

41.11 

10.84 

1.17 

0.01 

-48- 


PYRO  /  FILLED  CARTRIDGE  CASE 


INPUT  DATA 


PHI 

TEMP 

1  /T 

PHI/TSQUARE 

80.00 

504.67 

0. 19S1E-02 

0.31 41E-G3 

40  oOO 

495.67 

0.201 7E-02 

0. 1628E-03 

20.00 

489.00 

0 . 2045E-02 

0.8364E-04 

10.00 

482.04 

0.2075E-02 

0. 4304E-04 

5.00 

476.19 

0 . 2 l COE -02 

0. 2205  E -0  4 

2.50 

463.85 

0 ,2 l 33E-02 

0.1 137E-04 

E  (KCAL/MOLE)  0.445877E  02 
A  (SEC- 1 )  0.253103E  19 


DEGREES 

K 

PERCENT, 

UNREACTED,  (TIME  IN  YEARS) 

KELVIN 

(SEC  -1  ) 

1 .0 

2.0 

5.0 

10.0 

20.0 

298.0 

0.5  IE- 14 

100.00 

100.00 

100.00 

100.00 

100.00 

323.0 

0. 1 7E-1  1 

99,99 

99.99 

99.97 

99.95 

99.89 

348.0 

0.26E-09 

99.20 

98.40 

96.05 

92.25 

85.  1  1 

373.0 

0. 19E-07 

54.48 

29.68 

4.80 

0.23 

0.00 

PYRO  /  ULTRABOND  76-125 


INPUT  DATA 
PHI  TEMP 


80.00 

495.78 

40.00 

488.73 

40.00 

487.24 

20.00 

488.50 

20.00 

486.03 

10.00 

477.91 

10.00 

481.64 

5.00 

472.37 

5.00 

474.03 

2.50 

466.81 

1  /T 

0.20J7E-C2 
0.2046E-02 
0.2052E-02 
0.2047E-02 
0.2057E-02 
0 • 2092E -02 
0.2076E-02 
0.211 7E-02 
0.21 10E-02 
0,21 42E-02 


PHI /TSQUARE 
0.S255E-03 
0. 1  675E-03 
0. 1685E-03 
0.8381E-04 
0.3466E-04 
0.4378E-04 
0.431  IE-04 
0.2241 E-0  4 
0.2225E-04 
0.;  :  17E-C4 


E  (KCAL/MOLE)  0.526I31E  02 
A  (SEC-1 )  0. 199368E  23 


DEGREES 

K 

PERCENT, 

UNREACTED,  (TIME  IN  YEARS) 

KELVIN 

(SEC  -1) 

1.0 

2.0 

5.0 

10.0 

20.0 

293.0 

0.51E-16 

100.00 

100  ,00 

100.00 

100.00 

100.00 

323.0 

0.5CE-13 

100.00 

100.00 

100.00 

•00.00 

IOC. 00 

348.0 

0.18E-10 

99.94 

99.89 

99.72 

9S.43 

95.37 

373.0 

0.29E-08 

91.12 

83  .02 

62.30 

39.4 

15,56 

-52- 


i  •• 

i  !  ; 

;  I 

i 

.0019 

.0020 

.0021 

1/T 

~5  'i  - 

.0022 

![ 

_ _ I _ _ __ 

PYRO  /  EC  1099 


INPUT  DATA 


PHI 

TEMP 

1/T 

PHI/TSQUARE 

80.00 

501.97 

0. 1992E-02 

0.3175E-03 

40.00 

495.91 

0,201 6E-02 

0. 1627E-03 

40.00 

493.43 

0 . 2027  E-02 

0 . 1 643E-03 

40.00 

493.98 

0.2024E-02 

0. 1639E-03 

20.00 

487.43 

0.2052E-02 

0.841  8E-04 

20.00 

490.43 

0.2039E-C2 

0 . ZZ  1 5  E-0  * 

20.00 

488.30 

0.2048E-02 

0.8388E-C4 

10.00 

480.89 

0.2079E-02 

0.4324E-04 

10.00 

431.87 

0.2075E-02 

0. 4307E-04 

10.00 

481.84 

0.2075E-02 

0. 4307E-04 

5.00 

476.49 

0.2099E-02 

0 . 2202E-0  4 

5.00 

475.21 

0 . 2 1 0  4  E  -0  2 

0.221 4E-04 

5. CO 

474.53 

0 . 2 1 0  7  E  -0  2 

0 . 2220E-0  4 

2.50 

467.83 

0.2138E-02 

0.1 1 42E-0  4 

2.50 

463.76 

0.21 33E-02 

0.1  138E-04 

2.50 

468.32 

0.21 35E-02 

0.1  1  40E-04 

1.25 

461.74 

0 .2 1 66E-02 

0 . 5863E-05 

1.25 

463.68 

0 . 2 1  5  7  E  -0  2 

0.581 4E-05 

KCAL/MOLE) 

0. 465567E 

02 

SEC-1 ) 

0.223539E 

20 

DEGREES 

K 

PERCENT, 

UNREACTED,  (TIME 

IN  YEARS) 

KELVIN 

(SEC  -1) 

1  .0 

2.0 

5.0 

10.0 

20.0 

298.0 

0.I6E-1 4 

100.00 

100. oc 

100.00 

100.00 

100.00 

323,0 

0.72E-12 

100.00 

100  .00 

99.99 

99.98 

99.95 

348.0 

0. 13E-09 

99.59 

99.18 

97.95 

95,94 

92.05 

373.0 

C.12E-07 

68.56 

47.01 

15.15 

2.30 

0.05 

-■A- 


NACO  CONTROL 


’HI 

TEMP 

1/T 

PHI/TSQUi 

80.00 

501.27 

0. 1995E-02 

0.3 184E-03 

40.00 

492.44 

0.2031E-02 

0. 1650E-03 

40.00 

492.93 

0.2029E-02 

0 . 1 646E-03 

40.00 

491.70 

0.2034E-02 

0. 1654E-03 

20.00 

486.44 

0 . 2056E-02 

0 . 8452E-04 

20.00 

485.44 

0.2060E-02 

0 . 8487E-04 

20.00 

486.72 

0.2055E-02 

0 , 8443E-04 

10.00 

479.90 

0.2084E-02 

0.4342E-04 

10.00 

479.41 

0.2086E-02 

0.4351E-04 

10.00 

480.86 

0.2080E-02 

0 . 4325E-04 

5.00 

474.04 

0.21 IOE-02 

0 . 2225E-0  4 

5.00 

475.02 

0.210*E-02 

0.2216E-04 

5.00 

473.64 

0.21 1  IE-02 

0.2229E-04 

2.50 

466.37 

0.21 44E-02 

0.1 149E-04 

2.50 

467.64 

0.2138E-C2 

C.l 143E-04 

2.50 

467.39 

0.2140E-02 

0. 1 144E-04 

1.25 

460.29 

0.2173E-02 

0.5900E-05 

1.25 

462.23 

0.21 63E-02 

0.5851 E-05 

E  (KCAL/MOLE)  0.473649E  02 
A  (SEC-1)  0 . 643802E  20 


DEGREES 

K 

PERCENT, 

UNREACTED,  (TIME  IN  YEARS) 

KELVIN 

(SEC  -1) 

1.0 

2.0 

5.0 

10. 0 

20.0 

298.0 

0.12E-14 

100,00 

100.00 

100,00 

100.00 

100.00 

323.0 

0.57E-12 

100,00 

100.00 

99.99 

99.98 

99.96 

348.0 

0.1  IE-09 

99.64 

99.28 

98.20 

96.44 

93.01 

373.0 

0.1  IE-07 

69.96 

48.94 

16.76 

2.81 

0.08 

-56- 


NACQ  /  EC  1099 


INPUT  DATA 


PHI 

TEMP 

1/T 

PHI/TSQUARE 

80.00 

504.57 

0 . 1 9S2E-02 

0.31 42E-03 

40.00 

494.42 

0.2C23E-02 

0. 1636E-03 

40.00 

496.40 

0.2C15E-02 

0 . 1 623E-03 

40.00 

494.67 

0 

0. )  6351-03 

20.00 

483.45 

0.2068E-02 

0.8557E-04 

20.00 

485.44 

0.2060E-02 

0.84b7E-04 

20.00 

483.01 

0.2049E-02 

0.839SE-04 

10.00 

478.92 

0.2088E-02 

0.4360E-04 

10.00 

473.43 

0.2090E-02 

0.4369E-04 

10.00 

479.88 

0.2084E-02 

0.43  42E-0  4 

5.00 

473.55 

0.21 I2E-02 

0.2230E-04 

5.00 

475.51 

0.2103E-02 

C.221 IE-04 

5.00 

473.84 

0.21 10E-02 

0.2227E-04 

2.50 

467.34 

0.21 40E-02 

0.1 145E-04 

2.50 

466.86 

0.21 42E-02 

0.1 147E-04 

2.50 

466.71 

0.21 43E-02 

0.114SE-04 

1.25 

461 .26 

0.2168E-02 

0.5875E-05 

1.25 

461 .26 

0.2168E-02 

0 . 5875E-05 

£  (KCAL/MOLE)  0.437364E  02 
A  (SEC-1)  0.126858E  19 

DEGREES  K  PERCENT,  UNREACTED,  (TIME  IN  YEARS) 


KELVIN 

(SEC  -1  ) 

1.0 

2.0 

5.0 

10.0 

20.0 

298.0 

0. 1  IE-13 

100.00 

100.00 

100.00 

100.00 

100.00 

323.0 

0.32E-1  1 

99.99 

99.98 

99.95 

99.90 

99.80 

348.0 

0.43E-09 

98.65 

97.52 

93.44 

87.31 

76.23 

373.0 

0.30E  07 

39,01 

15.22 

0.90 

0.01 

0.00 

-58- 
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_7  _  ICtjAL  pel 
X  101*  per 

MOLE. 

SEC  OttO 

~l 

_ 

-3 

: 

STATISTICAL 
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_ 
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.  . 

_ 
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in 


59- 


NACO  /  ULTRABOND  76-125 


INPUT  DATA 


PHI 

ft  Mr 

1  /  i 

PHI/TSQUAR 

80.00 

501.97 

0. 1 992E-02 

0.3175E-03 

40.00 

491.94 

0.2033E-02 

0. 1653E-03 

40.00 

491 .94 

0.2033E-02 

0. 1653E-03 

40.00 

493.68 

0.2026E-02 

0. 1641E-03 

20.00 

484.94 

0 .2062E-02 

0.8505E-04 

20.00 

481.96 

C.2075E-02 

0.8610E-04 

20.00 

483.21 

0.2048E-02 

0.8391E-04 

10.00 

479.90 

0.2084E-02 

0 . 4342E-04 

10.00 

478.43 

0.2090E-02 

0 . 4369E-0  4 

10.00 

482.63 

0.2072E-C2 

0. 4293E-04 

5.00 

474.53 

0.2107E-02 

0 . 2220E-0  4 

5.00 

471.58 

0.2121E-02 

0.2248E-04 

5.00 

474.62 

C.2107E-02 

0 . 2220E-0  4 

2.50 

465.88 

0,21 46E-02 

0. 1 152E-04 

2.50 

465.39 

0.2149E-02 

0.1 154E-04 

2.50 

468.07 

0.21 36E-02 

0. 1 141E-04 

1.25 

459.32 

0,21 77E-02 

0.5925E-05 

1.25 

459.80 

0.2175E-02 

0.591 3E-05 

KCAL/MOLE)  0 

.  44629 7 E 

02 

SEC-1)  0 

.36 1038E 

19 

DEGREES  K  PERCENT,  UNREACTED,  (TIME  IN  YEARS) 


KELVIN 

(SEC  -1) 

1.0 

2.0 

5.0 

10.0 

20.0 

298.0 

0.67E-14 

100.00 

100.00 

100.00 

100.00 

100.00 

323.0 

0.23E-1 1 

99.99 

99.99 

99.96 

99.93 

99.36 

348.0 

0.34E-09 

98.94 

97.90 

94.83 

89.93 

80.87 

373.0 

0.25E-07 

44.81 

20,03 

1  .81 

0,03 

0.00 

-60- 


